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We compute the renormalization functions on the lattice, in the RI' scheme, of local bilinear quark 
operators \frF y/, where r = 1 , 75 , , 75 , 75 ay, v . This calculation is carried out to two loops for 
the first time. We consider both the flavor non-singlet and singlet operators. 
As a prerequisite for the above, we compute the quark field renormalization, Zy , up to two 

loops. We also compute the 1-loop renormalization functions for the gluon field, Z^ RI , ghost 

L RI^ L i?/' L, RI^ 

field, Z c ' , gauge parameter, Z a ' , and coupling constant Z g ' 

We use the clover action for fermions and the Wilson action for gluons. Our results are given as 
an explicit function of the coupling constant a = gl/16n 2 , the clover coefficient c$w> an d the 
number of fermion colors (N c ) and flavors (Nf), in the renormalized Feynman gauge. All 1-loop 
quantities are evaluated in an arbitrary gauge. 

Finally, we present our results in the MS scheme, for easier comparison with calculations in the 
continuum. We have generalized to fermionic fields in an arbitrary representation. Some special 
features of superficially divergent integrals, obtained from the evaluation of two-loop Feynman 
diagrams, are presented in detail in Ref. [[!]]. 
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1. Introduction 

Numerical simulations of QCD, formulated on the lattice, make use of a variety of composite 
operators, made out of quark fields. Matrix elements and correlation functions of a whole variety 
of such operators, are computed in order to study hadronic properties in this context. A proper 
renormalization of these operators is essential for the extraction of physical results from the lattice. 

In this work we study the renormalization function, Zp, of fermion bilinears & = i//Ty/ on 
the lattice, where r = i, 75, y^, YsYni Ys a nv (o]uv = 1/2[%>7v])- We consider both flavor singlet 
and nonsinglet operators. We employ the standard Wilson action for gluons and clover-improved 
Wilson fermions. The number of quark flavors Nf, the number of colors N c and the clover coeffi- 
cient csw are kept as free parameters. One necessary ingredient for the renormalization of fermion 
bilinears is the 2-loop quark field renormalization, Z v , calculated in The one-loop expression 
for the renormalization function Z g of the coupling constant is also necessary for expressing the 
results in terms of both the bare and the renormalized coupling constant. 

Our two-loop calculations have been performed in the bare and in the renormalized Feynman 
gauge. For the latter, we need the 1-loop renormalization functions Z a and Za of the gauge param- 
eter and gluon field respectively, as well as the one-loop expressions for Zp with an arbitrary value 
of the gauge parameter. 

The main results presented in this work are 2-loop bare Green's functions (amputated, one- 
particle irreducible (1PI)), for the scalar, pseudoscalar, vector, axial vector and tensor operator, 
as functions of the lattice spacing, a L , and the external momentum q. In general, one can use 
bare Green's functions to construct Z@ , the renormalization function for operator v, computed 
within a regularization X and renormalized in a scheme Y. We employ two widely used schemes 
to compute the various 2-loop renormalization functions: The RI' scheme and the MS scheme. 

The present work is the first two-loop computation of the renormalization of fermion bilinears 
on the lattice. One-loop computations of the same quantities exist for quite some time now (see, 
e.g., [J3j] and references therein). There have been made several attempts to estimate Zg non- 
perturbatively; recent results can be found in Refs. [Q-|6|]. A series of results have also been 
obtained using stochastic perturbation theory [f7|]. A related computation, regarding the fermion 
mass renormalization Z m with staggered fermions, can be found in [H]. 



2. Formulation of the problem 

We will make use of the Wilson formulation of the QCD action on the lattice, with the addition 
of the clover (SW) term for fermions. In standard notation, it reads: 



Sl = S G + ££(4r + m )v/ / (x)y / /-(x) 



1 



%(x) {r-Yv) U x , x +n Y f (x + n) + %(x + a0 {r + Yp) U x+fl , x yr f (x) 



(2.1) 



t c sw£ £ xj/f(x)a lxv F tlv (x)\if f (x), 
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1 , 

where : F^ v = ^ (Qfiv - Qv^i) (2-2) 
and: Q^y = U X}X +ixU x +p tX +p+vUx+p+v,x+vUx+v,x~\~U Xj x+vUx+v,x+v— nUx+v—n,x—nUx—[i,x 

+ U XjX —ji Ux-n^x-fx-v U x -fi,-v,x-v U x -v,x + U XjX — V U x -v,x-v+n U x -v+n,x+n ^x^\i,x (2.3) 

So is the standard pure gluon action, made out of 1 x 1 plaquettes. r is the Wilson parameter 
(set to r = 1 henceforth); / is a flavor index. Powers of the lattice spacing a L have been omitted 
and may be directly reinserted by dimensional counting. 

The "Lagrangian mass" m Q is a free parameter in principle. However, since we will be using 
mass independent renormalization schemes, all renormalization functions which we will be cal- 
culating, must be evaluated at vanishing renormalized mass, that is, when m is set equal to the 
critical value m cr : m -> m cr = mi g 2 Q + &(gi). 

One prerequisite to our programme consists of the renormalization functions, Z c , Z v , Z g 
and Z a , for the gluon, ghost and fermion fields (A^, c a , \\f), and for the coupling constant g and 
gauge parameter a, respectively (for definitions of these quantities, see Ref. [§]); we will also need 
the fermion mass counterterm m cr . These quantities are all needed to one loop, except for Zy which 
is required to two loops. The value of each Zg depends both on the regularization X and on the 

X Y 

renormalization scheme Y employed, and thus should properly be denoted as Z@ . 

As mentioned before, we employ the RI' renormalization scheme which is more immediate 
for a lattice regularized theory. It is defined by imposing a set of normalization conditions on matrix 
elements at a scale p., where (just as in the MS scheme): 

p. =n(4K/e^) l/1 (2.4) 

where y E is the Euler constant and jl is the scale entering the bare coupling constant g = jJ. e Z g g 
when regularizing in D = 4 — 2e dimensions. 



3. Renormalization of fermion bilinears 



The lattice operators Gy = must, in general, be renormalized in order to have finite 

matrix elements. We define renormalized operators by 



To 



(3.1) 



The renormalization functions Z r can be extracted through the corresponding bare 2-point 
functions L T (qa L ) (amputated, 1PI) on the lattice, through the employment of the RI 1 renormaliza- 
tion conditions: 



lim 

a L ^0 



7 L,RV 7 L,RV 

j v r 



(3.2) 



where L r (qa L ) is the appropriate bare 1PI 2-point Green's function on the lattice and r tree is its 
tree-level value. For the vector (V), axial-vector (AV) and tensor (T) operators, we can express the 
bare Green's functions in the following way: 



y^\qa L ) + q ^\qa h ) 
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^av(^l) = 757m4v(^) + 75^4?(<Z«l) (3-3) 

H 

MtfaJ = 75^v £ r W« L ) + 75 2 L t WhJ 



Only the 4av r (*7 a L ) P^ts are involved in Eq. ( |3.2| ). It is worth noting here that terms which break 



Lorentz invariance (but are compatible with hypercubic invariance), such as 7^ (q 11 ) 2 /q 2 , turn out 
to be absent from all bare Green's functions; thus, the latter have the same Lorentz structure as in 
the continuum. 

For easier comparison with calculations coming from the continuum, we need to express our 
ilts in the MS scheme. For each renormal 
its MS counterpart using conversion factors: 



results in the MS scheme. For each renormalization function on the lattice, Z@ RI , we can construct 



jL.0 gDRJtl' 

C r (g, a) = = r _ (3.4) 

1 V6 ' ; jLMS Z DR,MS V ' 

These conversion factors are regularization independent; thus they can be calculated more easily 
in dimensional (DR), rather than Lattice (L), regularization, (see, e.g., Ref. [[Kj). Due to the non- 
unique generalization of y 5 to D dimensions, the pseudoscalar and axial-vector bilinear operators 
require special attention in the MS scheme. 

For a more detailed analysis of the renormalization of fermion bilinears and their conversion 
to the MS scheme, see Refs. [fiL 0J. 

4. Computation and Results 

The Feynman diagrams contributing to the bare Green's functions, at 1- and 2-loop level, are 
shown in Figs. [T] and ||, respectively. For flavor singlet bilinears, there are 4 extra diagrams, shown 
in Fig. ||, which contain the operator insertion inside a closed fermion loop. These diagrams give a 
nonzero contribution only in the scalar and axial-vector cases. 



Figure 1: One-loop diagram conttibuting to Zp. A wavy (solid) line represents gluons (fermions). A cross 
denotes the Dirac matrices 1 (scalar), 75 (pseudoscalar), 7^ (vector), 757^ (axial vector) and 75 Cfyiv (tensor). 

In Figs. I tog' 'mirror" diagrams (those in which the direction of the external fermion line is 
reversed) should also be included. In most cases, these coincide trivially with the original diagrams; 
even in the remaining cases, they can be seen to give equal contribution, by invariance under charge 
conjugation. 

The evaluation of all Feynman diagrams leads directly to the corresponding bare Green's func- 

T L RI 1 

tions Ep. These, in turn, can be converted to the corresponding renormalization functions Z r ' , via 



Eq. (3.2). One-loop results for Z r ' are presented below in a generic gauge. The errors appearing 



in these expressions, result from an extrapolation to infinite lattice. 

3 ln{a 2 pL 2 ) - a Q - 16.9524103(1) 
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7L.RI' 



1 + 



Z^' RI ' = 1 + 



16k 2 



16k 2 



CF 



CF 



-7.7379159(3) c sw + 1.38038065(4) c| w 
3 ln(a^ 2 ) - a Q - 26.5954414(1) 
+2.248868528(3) c sw - 2.03601561 (4) c 2 sw 

-ln(a 2 jQ 2 ) + a - 17.018079209(7) 
+3.91333261(4)c sw + 1.972295300(5) c 2 sw 



(4-1) 



(4.2) 



(4.3) 



The corresponding expressions for Z v ' , Z AV can be read off from Eqs. (4.4,4.5) below. 

We present below Z v ' and Z A y to two loops in the renormalized Feynman gauge. The 
corresponding plots are exhibited in Figs. ^ and ||, as functions of the clover parameter, csw- 
For a complete set of results, regarding the renormalization functions and renormalized Green's 
functions, both in the RI' and in the MS scheme, the reader should refer to Refs. |[lj, g]. Furthermore, 
a calculation regarding the multiplicative mass renormalization Z m , which is directly related to 
the flavor singlet scalar operator, can be found in Ref. [^]. The generalization of our results to 
an arbitrary representation, as well as a detailed discussion regarding the superficially divergent 
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Figure 2: Two-loop diagrams contributing to Zp. Wavy (solid, dotted) lines represent gluons (fermions, 
ghosts). A solid box denotes a vertex from the measure part of the action; a solid circle is a mass counterterm; 
crosses denote the Dirac matrices 1 (scalar), 75 (pseudoscalar), 7^ (vector), 757^ (axial-vector) and Ys^nv 
(tensor). 




12 3 4 

Figure 3: Extra two-loop diagrams contributing to Z^_ s i n ^ et and Z^y singlet- A cross denotes an insertion of a 
flavor singlet operator. Wavy (solid) lines represent gluons (fermions). 
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integrals, can also be found in these papers. 



-yL.Rl' 



1 + 



167T 2 
(167T 2 ) 



-20.617798655(6) +4.745564682(3)c sw +0.543168028(5)c sw 
; c F \N f (25.610(3) - 11.058(1) c sw +33.937(3) ci w 



- 13.5286(6) c| w - 1.2914(6) <V 
+c F (-539.78(1) -223.57(2)c sw - 104.1 16(5) cl w 

-32.2623(8) cj w + 4.5575(3)4^ 
+N C ( - 5 1 .59(1) + 18.543(5) c sw + 20.960(6) c 2 sw 

+2.5121(5)c| w +0.1765(l)^ w 



(4.4) 




Figure 4: Zy RI (a L fi) — Zy MS (a L jx) versus csw (N c = 3, /i = l/a L , j8 = 6.0). Results up to 2 loops are 
shown for Nf = (solid line) and Nf — 2 (dashed line); one-loop results are plotted with a dotted line. 



7 L,RI' _ , 



167T 2 

.4 



CF 



15.796283066(5) - 0.247827627(3) c sw + 2.251366176(5) cj w 



+ ^2j2 c f [ N f (l8.497(3) - 1.285(l)c sw + 19.071(3)c| w 
+ 1 .0333(6) c 3 sw - 6.7549(6) 



+c F ( - 184.01(1) -389.86(l)c sw - 166.738(6)c| w 

+7.894(1) c^ w +4.3201 (3) c| w ) 
+A^ C ( -21.62(1) -33.652(5) c sw + 26.636(6) c| w 

+ 10.2186(5)cl w + 1.4893(l)c^ w N 



(4.5) 
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Figure 5: Z^f 7 (a L £i) versus csw (N c = 3, p. = l/a L , j3 = 6.0). Results up to 2 loops, for the flavor 
nonsinglet operator, are shown for Nf = 0, 2 (solid line, dashed line); 2-loop results for the flavor singlet 
operator, for Nf = 2, are plotted with a dash-dotted line; one-loop results are plotted with a dotted line. 
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